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We study the angular structure of polarization of light transmitted through a 
nematic liquid crystal (NLC) cell by theoretically analyzing the polarization state as 
a function of the incidence angles. For a uniformly aligned NLC cell, the 4x4 matrix 
formalism and the orthogonality relations are used to derive the exact expressions for 
the transmission and reflection matrices. The polarization resolved angular patterns 
in the two-dimensional projection plane are characterized in terms of the polarization 
singularities such as C points (points of circular polarization) and L lines (lines of 
linear polarization). In the case of linearly polarized plane waves incident on the 
homeotropically aligned cell, we present the results of detailed theoretical analysis 
describing the structure of the polarization singularities. We apply the theory to 
compute the polarization patterns for various orientational structures in the NLC 
cell and discuss the effects induced by director orientation and biaxiality. 
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I. INTRODUCTION 

The optical properties of nematic liquid crystals (NLCs) have long been known to play a 
key part in a large variety of their applications [TJ [21 Typically, NLCs are used as optically 
anisotropic materials where the anisotropy is determined by the orientational structure which 
is sensitive to external fields and, in restricted geometries, can also be influenced by changing 
the boundary conditions. 

Similar to other anisotropic materials, the polarization state of light propagating through 
a liquid crystal varies due to the presence of the anisotropy. When the NLC cell is placed 
between two crossed polarizers, these anisotropy induced changes of the polarization manifest 
themselves in variations of the intensity of the transmitted light. This effect underlines the 
vast majority of experimental techniques devised to characterize orientational structures in 
NLC cells. 

For example, the crystal rotation technique [I] relies on analyzing the angle dependence 
of the transmittance of NLC cells placed between crossed polarizers. It can be regarded as 
a special case dealing with one- dimensional cross-section of the conoscopic pattern. 

Nowadays the conoscopy is proved to be a useful tool for studying liquid crystal systems. 
It is used to detect biaxiality of NLCs [5l [6] and to measure the pretilt angle in uniaxial liquid 
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crystal cells [7J [8]. Orientational structures and helix unwinding process in ferroelectric 
smectic liquid crystals were also studied by conoscopy in Refs. [HI [10]. The conoscopic 
patterns of hybrid NLC cells were investigated in Ref . [Uj . 

In this paper we are aimed to study the polarization structure behind the conoscopic im- 
ages by performing a comprehensive analysis of the polarization state of the light transmitted 
through a NLC cell as a function of the incidence angles. Technically, we deal with angular 
two-dimensional distribution of the Stokes parameters describing the field of polarization 
ellipses which might be called the polarization resolved conoscopic (angular) pattern. 

It was originally recognized by Nye [121 DSl EH] that the important elements characterizing 
the geometric structure of such polarization fields are the so-called polarization singularities. 
In particular, these are the C points [the points where the light wave is circularly polarized] 
and the L lines [the curves along which the polarization is linear]. Over the past two 
decades these singularities and related issues have been the subject of numerous theoretical, 
experimental and numerical studies [15l [16l HZl UHl [19l [2QI EH E21 [23l EH [25l [26l E7j . 

The theory of polarization singularities was applied to study the angular dependence of 
the polarization of plane wave eigenmodes in birefringent dichroic chiral crystals [28]. The 
polarization state of the electric displacement field was analyzed in relation to the direction of 
the wave vector. This analysis was recently generalized and extended to a more complicated 
case of bianisotropic media [29] . 

The experimental results and theoretical analysis presented in Ref. [3U] deal with the 
unfolding of a linearly polarized Laguerre-Gauss (LGoi) beam with an on-axis vortex on 
propagation through a birefringent crystal. It was found that a complicated pattern of 
polarization singularities is formed as a result of the anisotropy induced symmetry breaking. 

Paraxial propagation of a beam carrying the angular momentum along the optic axis 
of a uniaxially anisotropic medium was also studied in Refs. [3U [321 E3]. Interestingly, 
there are a number of dynamical effects induced by a light beam normally impinging on a 
homeotropically aligned NLC cell [Ml ESI [361 EH EE] where liquid crystals appear as unique 
anisotropic materials sensitive to the spin and orbital angular momentum of light. 

In this study we adapt a systematic theoretical approach and explore the characteristic 
features of the polarization structure (the polarization resolved angular pattern) emerging 
due to interference of four eigenmodes excited in NLC cells by the plane waves with varying 
direction of incidence. 

The layout of the paper is as follows. The problem of light transmission through a 
uniformly anisotropic NLC cell is considered in Sec. [TT} Using the 4x4 matrix formalism 
and the orthogonality relations we deduce the exact expressions for the transmission and 
reflection matrices. In Sec. Ill the analytical results are used for analysis of the polarization 
resolved conoscopic patterns when the incident light is linearly polarized. 

The case of homeotropic NLC cells is treated in Sec. IIIB| We show that there are concen- 
tric circles each containing four symmetrically arranged C points of alternating handedness. 
For these C points, we obtain the analytical expressions for the topological index and for the 
discriminant whose sign determines the type of the point according to the line classification. 
The index and the handedness of the C points are found to alternate in sign along the radial 
direction. Morphology of the C points lying on the circle of the smallest radius representing 
the direction close to the normal incidence is shown to be of the lemon type with the index 
equal to +1/2. In sufficiently thick cells, other C points are either of the star type (the index 
is —1/2) or of the monstar type (the index is +1/2). The L lines are formed by two rotated 
coordinate axes and by circles separating the circles with the C points. Some numerical 
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results for the tilted, planar and biaxial structures are presented in Sec. HI C 
Discussion and concluding remarks are given in Sec. |IV| 



II. TRANSMISSION OF LIGHT THROUGH LIQUID CRYSTAL CELLS 

We consider a nematic liquid crystal (NLC) cell of thickness d sandwiched between two 
parallel plates that are normal to the z axis: z = and z = d. The NLC represents a 
birefringent material with the dielectric tensor given by 

e = e 3 Z 3 + Aeid<g>d + Ae 2 m <g> m, (1) 

where Ae, = 6j — e 3 and X n is the n x n identity matrix. (In what follows carets will denote 
unit vectors.) 

Typically, anisotropy of nematics is locally uniaxial and NLC molecules align on average 
along a local unit director [T]. In this case the NLC director d determines the optic axis 
and the expression taken in the limit of vanishing biaxiality with Ae 2 = gives the 
uniaxially anisotropic dielectric tensor. Its two principal values e 3 = e± and ei = e\\ define 
the ordinary and extraordinary refractive indices, n Q = y/fie± and n e = ^/yuey, where /i is 
the NLC magnetic permeability. 

In a more general case of biaxial nematics [391 HQ] that were recently observed experi- 
mentally [HI E], there are three different dielectric constants ei, e 2 and e 3 representing the 
eigenvalues of the dielectric tensor (II]), so that the eigenvectors d, m and 1 = d x m give the 
corresponding principal axes. The unit vectors d, m and 1 can be conveniently expressed in 
terms of Euler angles as follows 

d = sin 6*d cos 0d x + sin 6 A sin d y + cos #d z , (2a) 
m = cos 7d e x . (d) + sin 7d (d) , (2b) 
1 = -sin7 d e a: (d) + cos7 d e y (d), (2c) 

where e x (d) = (cos#d cos0d, cos^d sin0 d , — sin^d) and e y (d) = (— sin d , cos 4>di 0). 

We shall need to write the Maxwell equations for a harmonic electromagnetic wave (time- 
dependent factor is exp{— iut}) in the form: 

V x E = z/i£; vac H, (3a) 

V x H = -ik vac B, (3b) 

where k vac = u/cis the free-space wave number; /i is the magnetic permittivity and D = e-E 
is the electric displacement field. The medium surrounding the NLC cell is assumed to 
be optically isotropic and characterized by the dielectric constant e m and the magnetic 
permittivity /i m . So, Maxwell's equations in the medium outside the cell can be obtained 
from Eq. ^ by replacing [i and D with fx m and e m E, respectively. 

Referring to Fig. [TJ there are two plane waves in the lower half space z < bounded by 
the entrance face of the NLC cell: the incoming incident wave {E inc ,H inc } and the outgoing 
reflected wave {E re fi, H ro fi}. The transmitted plane wave {E tr , H tr } is excited by the incident 
wave and propagates along the direction of incidence in the upper half space z > d after the 
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FIG. 1: Geometry of nematic cell in the plane of incidence. 



exit face. So, the electric field outside the cell is a superposition of the plane wave solutions 
of the Maxwell equations 

EUo = E inc (k inc y( kl -- r ) + E refl (k refl )e i ( k -«- r ), (4a) 
E\ z>d = E tr (k tr )e l ( k - r \ (4b) 

where the wave vectors ki nc , k re fi and k tr that are constrained to lie in the plane of incidence 
due to the boundary conditions which require the tangential components of the electric and 
magnetic fields to be continuous at the boundary (interface) surfaces. These conditions are 
given by 

-P(z) • [E|, =0+0 - E| z=0 _ ] = V(z) ■ [E\ z=d+0 - E| 2=d _ ] = 0, (5a) 
P(z) • [H| z=0+0 - H| z=0 -o] = V(z) ■ [U\ z=d+0 - H\ z=d _ ] = 0, (5b) 

where "P(z) = 1% — z ® z is the projector onto the plane with the normal directed along the 
vector z (the x-y plane). 

Another consequence of the boundary conditions ^ is that the tangential components 
of the wave vectors are the same. Assuming that the incidence plane is the x-z plane we 
have 

k a = k m k a = k x x + z, a G {inc, refl, tr}, (6) 
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where k m /k v & c = n m = y/fj, m e m is the refractive index of the ambient medium and the 
components can be expressed in terms of the incidence angle 9- inc as follows 

k x = ^msin#i nc , (7) 

k (inc) = k (tr) = _ k (rcR) = ^ ^ ^ (g) 

The plane wave traveling in the isotropic ambient medium along the wave vector ^ is 
transverse, so that the polarization vector is given by 

E a (k a ) = E^e x {k a ) + E^e y (k a ), (9) 
e x (k a ) = k- 1 (k^SL-k x x) 1 e y (k a ) = y, (10) 

where and E^ are the in-plane and out-of-plane components of the electric field, 
respectively. The vector characterizing the magnetic field is 

/i m H a (k a ) = q a x E a (k a ) = n m [E^y - E { ° ] e x (k a )], (11) 
where q Q = k~^ c k a = n m k a . 



A. Operator of evolution and transmission matrix 



The expressions (J4|) — ( 1 1 ) give the electromagnetic field of incident, transmitted and re- 
flected waves propagating in the ambient medium. This field is of the general form 



E{x,z) = E(z)e 



H(x,z) = HUe ,Sl 



(12) 



that can also be used for the field inside the NLC cell when the dielectric tensor Q is 
independent of the in-plane coordinates x and y. 

The boundary conditions ^ can be conveniently written using matrix notations as follows 



where 



F< 
F> 



O | Einc 



E 




tr 



= F(0), 
F(d), 



(13a) 
(13b) 



F(z) = 



( E x (z) \ 

Ey( Z ) 



E a = 



E \ 
E 



(14) 



Our notations can be regarded as a version of the well known 4x4 matrix formalism due 
to Berreman |42j . 

The matrix Q m relates the field vectors, F< = F(0 — 0) and F> = F(d + 0), and the 
vector amplitudes E a of the waves in the surrounding medium. Its block structure 



7i T-C^ 



(15) 
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is characterized by the four 2x2 matrices: and "H.± . From Eqs. (|9j)-( 11 ) it can be 

readily obtained that these matrices are diagonal and are given by 



/^m 



(16a) 
(16b) 



where cr 3 = diag(l, — 1), A = diag(l, cos# inc ) and B = /W^n 1 diag(cos 9i nc , 1). It is rather 
straightforward to check the validity of the orthogonality relations for the matrices (16) 



(17) 



where J\f+ 



±2^ TO 1 n m cos0i nc X 2 , is the Kronecker symbol and the superscript f 



indicates Hermitian conjugation. 

As it was shown in Ref. [13], the relations (17), where A/" m is a diagonal matrix, can be 
deduced from the conservation law for the energy flux in non-absorbing media that, in the 
case of one-dimensional inhomogeneity, requires the z component of the Poynting vector to 
be independent of z. The relations (17) can be rewritten in the more compact matrix form: 

A/" m , Q e 



[QjQQr 



o x 2 
x 2 o 



(m) K Km)\ 



iV m diag(X 2 , -X 2 ) and iV n 



where A/" m = diag(A/+" ,JV 

braically, Eq. (18) can be used to simplify inversion of the matrix Q n 

Q" 1 = (2 cosine)" 1 diag(X 2 ,cr 



A B 

-A B 



(18) 

?7. m cos^i nc . Alge- 
(19) 



and to ease qualitative analysis [see , e.g., Ref. jH] for a more extended discussion of appli- 
cations]. 

Given the characteristics of the incident wave, the amplitudes of the reflected and trans- 
mitted waves cannot be computed from the boundary conditions ( 13 ) unless the electromag- 
netic field within the NLC cell is defined. At this stage we just note that, for the in-plane 
components of the electric and magnetic field, {E x ,E y } and {H x ,H y }, the result can be 
presented in the form of a general solution to the linear problem 



F(z)=W(*)F(0), 



(20) 



where the matrix- valued function U{£) is the evolution operator which can be computed by 
solving the Maxwell equations ([3]). 

After substituting Eq. (20 ) into the boundary conditions (13a) and using the orthogonality 
relations (17), we have 



V 







v 



E re fl 



where U 1 = U~ l (d) and Q 



B -B 
A A 



Vhe Vhh 



diae(X9,cr ? 



(21) 
(22) 
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From Eq. (21) we have 



E l 
E 



(tr) N 
(tr) 



r 



E i 

E 



(mc) N 
II 

(inc) 



r = v 



-i 



(23) 
(24) 



where T is the transmission (transfer) matrix linking the transmitted and incident waves. 
Similar result for the reflected wave is 



E, 



(refl) N 



E 



(inc) x 



= v h y-} = v he r, 



(25) 
(26) 



where 1Z is the reflection matrix. 



B. Eigenmodes 



The equations (24) and (26) give the transmission and reflection matrices expressed in 



terms of the matrix (22) which is related to the evolution operator (20). For our purposes 



it is instructive to see how this operator can be derived by using the basis of eigenmodes 
(normal modes) in a uniformly anisotropic medium. In this case the eigenmodes are known 
to be linearly polarized plane waves characterized by the k dependent refractive indices [4"5l 



For plane waves, owing to the Maxwell equation (3b), the electric displacement field is 
transverse, (D • k) =0. So, assuming that the wave vector k lies in the plane of incidence 
(the x-z plane), the vector D can be conveniently defined by its components in the basis 
{e x (k),e y (k),e z (k) = k} 



D = D x e x {k) + D y e y {k), 

k = e 2 (k) = k~ x k = k~ (k x Sc + k z z) = g _1 q, 



(27) 
(28) 



where e x (k) = q 1 (q z x-q x z), e y (k) = y, q x>x = k XtZ /k vac and q = k/k VSiC = n is the refractive 



index. The vector q defined in Eq. (28 ) is parallel to k and its length is equal to the refractive 



index. For convenience, we shall often use this vector in place of the wave vector. 



Given the electric displacement field (27), the electric field can be found from the consti- 
tutive relation 



E = fit] ■ D, 

77 = Tjs X 3 + Ar)i d ® d + Ar/ 2 rh <S> rh, 
where rji = (/xej) -1 , Arji — r\i — 773 and \ir) is the inverse dielectric tensor [firj 



(29) 
(30) 



The Maxwell equations can now be combined with the relation (29) to yield the equation 



for the displacement vector D in the form of an eigenvalue problem 



r7, D = g- 2 D, m = V(k) ■ 77 ■ V(k), 



(31) 
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k <S> k is the projector onto the plane normal to the wave vector k. By 



where 'P(k) = X 

using the expression for the inverse dielectric tensor (30) we can rewrite the equation (31) 
in the explicit matrix form: 



[a_ cr 3 + bcr 1 



«•] ( D. 



0, (Tl 



1 

1 



2a T = Ar)i(d 2 x T d\) + A V2 {m 2 x =F m% 
b = Ar)id x d y + Arj 2 m x m y , 



A 



VaQ 



(32) 

(33) 
(34) 
(35) 



where d x<y = q(d- e X) y(lk)) and rh X) y = q(m-e Xty (k)). Then the dispersion relation (Fresnel's 
equation) 

(l- V3 q 2 )[l- V3 q 2 -A Vl (q 2 -d 2 z )- 
Ar] 2 (q 2 - m 2 z )] + Ar) X Ar) 2 q 2 l 2 z = 0, 



where ql z 

equations (|32|) to have a non-vanishing solution. 



(36) 

q 2 (l • k) , can be derived as the condition for the system of linear 



The Fresnel equation describes the wave surface. In our case solving the algebraic equa- 
sin# inc gives the values of the z component of the vector q, q z . 



tion (36) at q x 



Generally, there are four roots of Eq. (36) qi a \ Each root corresponds to the eigenwave 
propagating inside the cell with the wave vector k a = k vac (q x , 0, q^) and the refractive index 
n a = got- The corresponding polarization vector of the electric displacement field is given by 



where 



D («) _ j c os<p a e x {k a ) + sm(p a e y {k a ), X a > 0, 
\-sin(f) a e x (k a ) + cos0 a e y (k a ), X a < 0, 

^« = M gz=g (^i 2 Q = arg(a_ + i&) | „_ = „(«) ■ 



1z=q z 



(37) 



(38) 



From Eq. (38 ) it is clear that the azimuthal angle <p a becomes indeterminate in the degenerate 
case when the coefficients a_ and b are both identically equal to zero. Typically, as far as the 
eigenmodes are concerned, this case does not present any fundamental difficulties. It just 
means that the azimuthal angles of the degenerate eigenmodes can be prescribed arbitrarily. 
Such freedom of choice, however, does not affect the evolution operator which remains 
uniquely defined. 

The procedure to determine the characteristics of the eigenmodes involves the following 



steps: (a) evaluation of q z by solving the Fresnel equation (36 ); (b) calculation of the polariza- 
tion vectors of the electric displacement field by using the formula (37); (c) computing 



the polarization vectors of the electric and magnetic fields from the relations: E^) = /i Ty-D^ 
(see Eq. d29b) and /iH^ = q a x (see Eq. Q). 



1. Uniaxial anisotropy with Ae 2 = 



Now we apply the above procedure to the limiting case of uniaxial anisotropy with e 2 = e 3 
and Ae 2 = Ar] 2 = 0. At Arj 2 = 0, the Fresnel equation (36) takes the factorized form and 
the values of q z can be found as roots of two quadratic equations. 
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The first equation 1 — rj^q 2 = represents the spherical wave surface. The corresponding 
eigenmodes are known as the ordinary waves. There are two values of q z 



Q. 



r ) = ±v»i-& (39) 



where n\ = [163 = fie±, that are equal in value but opposite in sign. When, similar to the 
incident and transmitted waves, the z component of the wave vector (and the vector q) is 
positive, the eigenmode might be called the refracted (forward) eigenwave. In the opposite 
case where, similar to the reflected wave, qz is negative, the eigenmode will be referred to 



as the reflected (backward) eigenwave. So, Eq. (39) describes two ordinary eigenmodes: the 
refracted eigenwave with q z = qi + °^ > and the reflected eigenwave with q z = q z ~°^ < 0. 
The second equation 

q 2 + u a (q-d) 2 -nl = 0, (40) 

where njj = /iei = //ey and u a = —Arji/rjx = (n| — n\)/n\ is the anisotropy parameter, gives 
the values of q z for the eigenmodes known as the extraordinary waves. These are given by 



qf e) = [1 + Uad 2 ]- 1 {-u a d z d x q x ± ^D}, (41) 
D = nj (1 + u a d 2 z ) -q 2 x [l + u a (d 2 x + dl)] (42) 



where d x = (d • x) and d z = (d • z) . 



At u a > (u a < 0), in the x-z plane, Eq. (40) describes the ellipse with the major (minor 



semi-axis of the length n\\ oriented perpendicular to the projection of the director (2a) on 



the plane of incidence (d x ,0,d z ). The length of minor (major) semi-axis, n± = [n~ ± 
u a{d y /n\\) 2 ]~ l l 2 , depends on the y component of the director and varies from n± to n\\ as d 2 
increases from zero to unity. Clearly, degeneracy in refractive indices with n Q = n± e may 
occur only if the director is in the incidence plane (0d = 0). Additionally, the matching 
condition for the x components of q and the director q x = n m sin ^ inc = ±n Q d x = ±n sin 
needs to be met. 

The wave vectors and the refractive indices of the normal modes are determined by the 
relation 

q± a = ^~ac k ±a = Qx x + q[ ±a) z = n ±Q k ±a , a e {o, e}, (43) 

where n± a = q± a , n± = n = n± is the ordinary refractive index and n± e is the refractive 
index of the extraordinary wave propagating along the unit vector k_|- e . 
The relations 

sign A io = - sign A771 , sign A ±e = sign A771 , (44) 
a_ + ib = Arji(d x + id y ) 2 (45) 



obtained by substituting A7/2 = into Eqs. (33)-(35) can now be combined with Eq. (37) 



to yield the following result for the polarization vectors of the electric displacement field: 

D (±o) oc -d y e x (k± e ) + d x e y (k± e ) oc k io x d, (46) 
oc d x e x (k ±e ) + d y e y (k ±e ) oc V(k ±e ) ■ d. (47) 
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Following the procedure described at the end of Sec. II B we then find the polarization 
vectors of the electric field for the eigenmodes 



E (±o) = -q io x d, 

E (±e) = [d-^ 2 (d-q ±e )q ±e ], 



where Eq. ( 48b ) was derived by using the identity 

77 ■ P(k)d = n- 2 d - nl 2 (d ■ k)k 



(48a) 
(48b) 

(49) 



.2] _ 2 



with n 2 = q 2 = [l + u a (d • k) ] ; 

The result for the magnetic field of the normal modes is 

H ( ± °)= / i- 1 [^d-(d-q ±0 )q ±0 ]. 
H (±e) = /i' 1 q ±e x d. 



(50a) 
(50b) 



In order to present the results in the matrix form, we, by analogy with Eq. (15), define 
the matrix composed of the polarization vectors as follows 



where 



Q 



E (±o) E (± e y 



E 



(±0) E (±e) \ > 



£+ S 



(51) 



(±0) 



H, 



(±e) 



rr(±o) „(±e) 
-Six —ti x 



(52) 



The explicit expressions for the in-plane components of the electric field are given by 



4 ±e) = n~ 2 {(n 2 - n 2 ±e )d x + # e) [q ±e x d]J, = d y . 

Those for the magnetic field are 

/itf ( ±0 > = q^ [q io X d] y , pLHW = n 2 d y , 
fiH^ = -qi ±e) d y} = [q ±e X d] y , 

where [q a xd] 9 = qi a) d x - q x d z . 



(53a) 
(53b) 

(54a) 
(54b) 



The matrix of eigenvectors (51) along with the eigenvalues given in Eq. (39) and Eq. (41) 
can now be used to describe the z dependence of the in-plane components of the electro- 
magnetic field inside the cell in terms of the eigenmodes amplitudes. So, we have 



F(z) = QUAz) 



0n 
/3- 



U d {z) 





(-)/ 



Uf\z 



where /3 d 



(3±o 



ur j wy (55) 

di&g(exp[iki ±o) z} J exp[iki ±e) z}) J (56) 
are the columns representing the amplitudes of the eigenmodes. 
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2. Transmission and reflection matrices 



Now it is not difficult to derive the evolution operator (20) from the equation (55). We 
can also use the orthogonality relation (18) for the matrix (51 ) so as to obtain the following 
result 



U(z) = QU d (z) Q 1 = QU d {z)M- l Q ] g, 



(57) 



where the matrix Q is defined in Eq. (18). 

Owing to the orthogonality relations (17), the matrix Af = Q}QQ, = diag(A/+, A/*_) is 
diagonal and its non-vanishing elements 



(58) 



are proportional to the normal components of the Poynting vector of the eigenmodes. Note 
that the evolution operator is not unitary. It is clear from the identity 

U-\z) = W{z) = GU\z)Q, (59) 

where an asterisk indicates complex conjugation, that immediately follows from the expres- 
sion (57). 



The operator (57) can be substituted into the expression for the matrix V (22) to yield 
the transmission and reflection matrices given by Eq. (24) and Eq. (26), respectively. By 

for the matrix V 

V* = SsV+Ss, Qz = diag(X 2 , -X 2 ). 



using Eq. ( 59 ) , we have the relation for the matrix V 

V 1 



(60) 



The conservation law T'T + IZ'IZ = X 2 is a consequence of the identity (60) and we also 
arrive at the conclusion that the transmission matrix (24) is symmetric. 



The explicit formulas for the transmission and reflection matrices given by Eq. (24) and 
Eq. (26). can be obtained by substituting Eq. (57) into Eq. (22). After some straightforward 
algebraic manipulations we have 

T = 2/i m n" 1 cos 9 inc r" 1 , r = r+ - 
B + U + \Aj + BJU\Aj 



where 



n = a, 



T± = A ± U ± [A ± ]\ U ± = 
A± = BU ± + A£±, B± 



T-, 
-1 



(61) 
(62) 



= BTL± 



AS . 



(63) 
(64) 



Substituting the expressions for the electric and magnetic fields of the eigenmodes [see 
Eq. (48) and Eq. (50), respectively] into Eq. (58) gives the diagonal elements of the matrices 
Af+ and JV_ 



/^ ±0 = 2g( ±0 )[^-(q ±0 .d) 2 ]. 



±e 



2n~ 



d z (q± e • d) 



n 



±e 



^) + # e) K-(q ±e -d) 2 ) 



(65) 
(66) 
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The expressions (61)-(66) can now be used to write explicitly the elements of the trans- 



mission matrix. The matrices £± and 7i± are given by Eqs. (52)-(p4ft and their diagonal 



(non-diagonal) elements are zero provided the director is normal (parallel) to the plane of 



incidence. In these cases, where either d y 



1 or cL = 0, the matrix r and the transmission 



matrix are both diagonal 



diag(r x , r y ), T = diag(t x , t 



in- 



where t n 



I, 1 cos 9 



The simplest case occurs when the cell is homeotropically aligned and d z 
homeotropic director structure, the elements of the matrix r are given by 

t x =2/x m n" 1 cos 9 inc cos 5 e - sin 5 e x 

[(/im^o^m 1 cos6» inc ) 2 + (n" 1 ^) 2 ] , 

T y =2/i m n" 1 cos ^ cos S - %i o) ] -1 sin 5 Q x 
[cos 2 inc + (/i^g'"') 2 ] , 

where gi o) = q%, q^ = n o n |j" 1 y n y — 9s anci ^ = qi^Ka^d. 



(67) 
1. For the 

(68a) 
(68b) 



III. ANGULAR STRUCTURE OF POLARIZATION ELLIPSE FIELD 

In the previous section we deduced the expression for the transmission matrix of the NLC 
cell and now we pass on to discussing how the polarization properties of the transmitted 
light depend on the direction of the incident wave. 

This direction is specified by two angles: the incidence angle 9 inc and the azimuthal angle 
of the plane of incidence inc . Clearly, we need to replace the director azimuthal angle 
d with d — inc so as to have the Euler angles describing orientation of the director Q 
with respect to the incidence plane. So, the transmission matrix T[9 mC} 0d) obtained in 
Sec. [Tf] is changed to T(9 inc , 9^, 4>d — <j>mc)- Dependence of the polarization parameters of 
transmitted waves on the angles 9- mc and 4>- mc will be of our primary concern. 

We begin with introducing necessary parameters and notations. Much of this material 
can be found in standard textbooks such as [121 SJJ (see, e.g., Ref. [27] for a more recent 
and extended discussion). 

The polarization state of a plane wave can be conveniently characterized by using the 
components of its vector amplitude in the spherical basis as follows 

E(k) = £+e+(k) + £_e_(k) = e^°{p+ + tp^}, (69a) 
E± = 2~ 1/2 (£|| TiE±) = \E±\ exp{i0 ± }, (69b) 
p ± = 2- 1 / 2 {|£ + |±|£_|}e : : iJ/ (k), o = (0_ + + )/2, (69c) 

where \/ze±(k) = e^k) ± ie y (h); the angle 0o is variously known as the phase of the 
vibration [12] or the rectifying phase |27J; p + and p_ are the major and minor semiaxes 
of the polarization ellipse with the principal axes directed along the unit vectors e^(k) = 
cos(f) p e x {k) + smcj) p ey(k) and e^(k) = — sm(p p e x (k) + cos <p p e y (K), respectively. So, the 
angle 

P = (0_ - + )/2 = 2- 1 arg([E + ]*E_) = 2" 1 arg (£_/£+) (70) 
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specifies the orientation of the polarization ellipse and will be referred to as the azimuthal 
angle of polarization or the polarization azimuth. 

Another important characteristics of the polarization ellipse describing its eccentricity is 
the signed ellipticity parameter 



6 ell = | F | , | F r (71) 



\ E - 




E+\ 


\E_ 


+ 


E + \ 



This parameter will be referred to as the ellipticity and its sign defines the handedness of 
the ellipse. 

According to Ref. [15], the ellipse is considered to be right-handed (RH) if its helicity is 
negative, (p + x p_ ■ k) < 0, so that \E + \ < |£L| and e c n > 0. For the left-handed (LH) 
ellipse, (p + x p_ • k) > and e e n < 0. 

Experimentally, the characteristics of the polarization ellipse can be obtained by measur- 
ing the Stokes parameters related to the coherence matrix with the elements Ai a f3 = E a E/3* , 
where a, (3 € {||, _L}. In circular basis, this matrix can be written as a linear combination of 
the Pauli matrices. The coefficients of the combination are the Stokes parameters 

Mc = CM& = ( ]f fffi) = 2- 1 J2 ^ *u (72) 



=o 



where C = 2 x l 2 ^ j }j , cr = X 2 and cr 2 = ^ . Since the determinant of the 

coherence matrix vanishes, det Al = 0, the Stokes parameters lie on the four-dimensional 
cone Sq = Yli=i $h an d can parameterized as follows 

S = \E + \ 2 + \E_\ 2 = \E {1 \ 2 + \E ± \ 2 , (73a) 



^11 

1 2 



51 = 2 Re E* + E_ = \E {1 \ 2 - \E ± \ 2 = S o cos2x p cos20 p , (73b) 

5 2 = 2 Im E\E- = 2 Re E L E\ = S cos 2 Xp sin 20 p , (73c) 

5 3 = \E + \ 2 - \E_\ 2 = 2lmE L El = S sm2 Xp , (73d) 



where < (p p < n is the polarization azimuth (70) and — tc /A < Xp ^ 7r/4 is the elliptic- 
ity angle. Then, the relations expressing the ellipse characteristics in terms of the Stokes 
parameters are 

p = 2" 1 argS, 5 = ^ + ^2, (74) 
e c ii = -tan x P , X P = 2" 1 arcsin(S , 3 /S , ). (75) 



Similar to Refs. [271118], we have used Eq. (72) and Eq. (73d) to define the Stokes parameter 



S3 which is opposite in sign to that given in the book [45J. 

The important special case occurs when the wave is circularly polarized and \E U \ = 0, 
so that the phases <p v and <p p are indeterminate. This is an example of the polarization 



singularity that, according to Eq. (74), can be regarded as the phase singularity of the 
complex Stokes field, S = Si + iS 2 - The point where the polarization is circular and \E„\ = 
will be referred to as the C v point with e e u = v. In our case such points are characterized 
by the incidence angles at which the transmitted wave is circular polarized and E v =0. 

The case of linearly polarized wave with \E + \ = \E_\ provides another example of the 
polarization singularity where the handedness is undefined. The curves along which the 
polarization is linear are called the L lines. 
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A. Angular patterns: C points and L lines 



Similar to Eq. (23), the transmission matrix in circular basis 

T ,. = i ' \ =CT& 




(76) 



relates the circular components of the incident and transmitted waves, {£f c) ,£i inc) } and 



Since the transmission matrix (|6l|) is symmetric, the diagonal elements of 



the matrix (76) are equal, t ++ = t = (t xx + t yy )/2, whereas the non-diagonal elements 

(t xx — t yy )/2 =F t xy differ provided T is not diagonal. 



t 



The transmission matrix 



T(p, 0) = exp(-i(f) <r 3 )T c (p, 0) exp(z0 cr 3 ) 
p = rtan6> inc , = inc 



(77) 
(78) 



describes the conoscopic patterns on the transverse plane of projection, where p and <ft 
are the polar coordinates (the Cartesian coordinates are: x = pcos0 and y = psin0) 
and r is the aperture dependent scale factor. We concentrate on the case, in which the 
incident plane wave is linearly polarized along the unit vector: cos ip p e ;c (ki nc )-|-sin ip p e y (k inc ). 



For the incident wave with the circular components E, 
components of the transmitted wave are given by 



(inc) 
v 



exp(-z^ P 



| S inc |, the reduced 



El tr) /\E inc \ = q u /2 = [t v , v + t Vt -„exp{-2ivif>j\ exp(-iz/^ ; 



vi- 



where ip = <fi — ip p . 

The polar coordinates of the C„ points on the projection plane, p£ and 
the numbering label, can be found by solving the equation 

l^(p,0)l = O 

that generally has multiple solutions. 

The C v points can be viewed as the phase singularities of the complex scalar field 

S = =§! + i§ 2 



(79) 



where k is 



(80) 



(81) 



proportional to the Stokes field defined in Eq. (74). Such singularities are characterized 



by the winding number which is the signed number of rotations of the two-component field 
(Si, S2) around the circuit surrounding the singularity [30]. The winding number also known 
as the signed strength of the dislocation is generically ±1. 

Since the polarization azimuth ( 74 ) is defined modulo 7r and 2<f) p = arg S, the dislocation 
strength is twice the index of the corresponding C u point, Ic- For generic C points, Ic = 
±1/2 and the index can be computed from the formula 



sign 



lm(d x S*d y S) 



(82) 



where d x f is the partial derivative of / with respect to x. 



The relation (82) gives the index of the C u point with the coordinates (x u ,y u ) expressed 



in terms of the vorticity [20], |2T] : lm(d x S*d y S) = d x Sid y S 2 — dySid x S 2 - The formula linking 
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gradients of the complex field (E ■ E) oc and the index for C lines in the three- 

dimensional space was derived in Ref. [50] . 

If ^ u = (C„ point), only derivatives of enter the expression (82) which can be 

(83) 



suitably rearranged to yield the index of the C u point in the following form: 



sign 



v 



sign 



Im(c^c^:) 



P=Pv 



Subsequently, we shall apply the formula (83) expressing the index in terms of the derivatives 



with respect to polar coordinates to the case of homeotropically aligned cell. 

In addition to the handedness and the index, the C points are classified according to the 
number of straight lines terminating on the singularity. This is the so-called line classification 
that was initially studied in the context of umbilic points [5T] ■ 

For generic C points, the number of the straight lines, Nc, may either be 1 or 3. This 
number is 3 provided the index equals —1/2, Ic = —1/2, and such C points are called stars. 
At Ic = 1/2, there are two characteristic patterns of polarization ellipses around a C point: 
(a) lemon with Nc = 1 and (b) monstar with Nc = 3 [12] . The quantitative criterion to 
distinguish between the C points of the lemon and the monstar types was deduced in [27J. 

We conclude this subsection with the remark that the transmitted wave is linearly polar- 
ized when the condition 



|vi/ + (p,0)I = I*-(p,0)I 



(84) 



is satisfied. So, Eq. (84) describes loci of points forming the L lines lying in the projection 
plane. 



B. Homeotropic cell 



When the director is normal to the substrates, the NLC cell is homeotropically aligned 
and d z = 1. In this case the transmission matrix is diagonal and the circular components of 
the transmitted light are given by 



[(t X + ty) + (i, 



ty) exp(— 2ii/^)] exp(— ivip 



pi- 



(85) 



where the transmission coefficients t 
angle of the incidence plane <fi = 



t 



x,y\ 



tx, y (p) do not depend on the azimuthal 



and are defined by Eq. (67) and Eq. (68). 



From Eq. (85), the C points may appear at the incidence angles 9^ 
condition: 



\tx ~t~ tn 



\t X ty\ 



Another form of this condition 



R(p) = Re(t x t*) = \t x \ ■ \t y \ cos 5 = 0, 



that satisfy the 



(86) 



where 5 is the phase difference, defines the radii of circles containing the C points. Note that, 
in weakly anisotropic NLCs, the phase difference 6 can be approximated by the phase shift 
due to the difference in optical path of the ordinary and extraordinary waves, 5 ~ 5 e — 5 , 
for the incidence angles up to 70 degrees. Clearly, Eq. (86) implies that 5 = 7r/2 + 7r/c, where 
k — 0, 1 ... Af — 1 is the non-negative integer and iV is the number of solutions. So, in the 
transverse projection plane, there are iV circles with C points. The radii of the circles pk, 
k 



0,1 ... N — 1, can be found by solving Eq. (86) 
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(c) Tilted structure at 9^ = 60 deg 



(d) Planar structure 



FIG. 2: The angular structure of polarization computed as the polarization ellipse field in the 
projection plane where the polar coordinates are p = tan#i nc and <f> = 0i nc . The incident light of 
the wavelength 633 nm is linearly polarized with ip p = 0. Other parameters used in calculations are: 
the cell thickness is d = 18.7 fim; \x = /i m = 1; n m = 1.5; n e = n\ = 1.71 and n Q = ni = 713 = 1.527 
(5CB). The C points of the lemon and the monstar types are marked by diamonds and triangles, 
respectively. Stars are used for the C points of the star type. Left-handed and right-handed 
polarization is respectively indicated by open and filled ellipses. 



At p = pfc, it is not difficult to obtain the expressions for the amplitude of the compo- 



nents (85) 



IMP* 



\t x \ z + \t 



2\t\ 2 [l + cos2{^ - v(-l) k a}]. 



tan a 



\tx 



(87) 
(88) 



From Eq. (87), it immediately follows that there are two pairs of C points with the azimuthal 
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(b) Biaxial structure at 0<j = 10 deg. 



FIG. 3: The angular structure of polarization computed as the polarization ellipse field in the 
projection plane where the polar coordinates are p = tan#i nc and <f> = 4>- mc . The incident light of 
the wavelength 633 nm is linearly polarized with ip p = 0. Other parameters used in calculations are: 
the cell thickness is d = 18.7 ^m; p, =1; 7<j = 45 degrees; n\ = 1.71, ni = 1.52 and 713 = 1.53. The 
C points of the lemon and the monstar types are marked by diamonds and triangles, respectively. 
Stars are used for the C points of the star type. Left-handed and right-handed polarization is 
respectively indicated by open and filled ellipses. 



angles given by 



<j>f k = ± tt/2 + v(-l) k a 



(89) 



in each circle of the radius pk- The symmetric arrangement of the C points [see also Fig. 2(a) )] 
is a consequence of the symmetry relations 



\*v(p,il>)\ = |^(p,7r + V)|, \V V M)\ = |tt-„(p,7r-V0| 



(90) 



for the amplitudes (85). 



We can now substitute Eq. (85) into the expression for the index (83) to derive the result 



I c = --d p R(p) 



-I] 



p=Pk 



(91) 



where the second equality follows because R is an oscillating function of p and R(0) > 0. 

Interestingly, Eq. (91) relates the index of C points and derivatives of the transmission 
coefficients with respect to the incidence angle (d p f = cos 2 9^ c df/d9i nc ). The index is thus 
determined by the circle number k and alternates in sign starting from Iq = +1/2. So, the 
C points in the vicinity of the origin at p = po can either be lemons or monstars. 

Now it is our task to formulate the criterion to distinguish between these two types. By 



using Eq. (85) we find the complex Stokes field (81) and the polarization azimuth (74) 



S = 4e z ^P, P = \t x \ 2 cos 2 ip - \t y \ 2 sin 2 iff - iR sin2^, 



4e 2i *P, 
b p = 20 + arg P 



(92) 
(93) 
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for the homeotropic cell. 

Our next step is to consider the first order term of the power series expansion of P in 
Ax = Arcos0 r and Ay = Arsin0 r . This is the linear part of the expansion that can be 
conveniently expressed in polar coordinates as follows 



pAP « p [Ax d x P + Ay d y P] 
Ar [cos(0 r — 4>)pd p P + sin(0 r 



d*P]. 



(94) 



At a C point, the expression (94) represents the complex field P expanded to first order. 
The identity 



Re P = \t\ 2 cos(ip — a) cos(ip + a). 
eases computing the derivatives at the C point. The result is 



d^P 
pd P P 



P=Pk 
i/)=7r/2±a 



P=Pk 
ip=ir /2±q 



±|t| 2 sin(2a), 
= sin(2a) [— \t\ 2 cei ± iR\\ . 



Pi = {pd p R) 



p=Pk 



ai = (pdpa) 



p=pk 



(95) 

(96) 
(97) 
(98) 



We can now follow the line of reasoning presented in Ref. [27] and derive the polynomial 
equation 



tan 2((f) p — 



2q 



Ri 



l- g 2 \t\2( q if ai ) 



ImAP 
Re AP 



(99) 



in q = tan(0 p — 0) = tan(0 r — 0). The roots of Eq. (99) give the angles where P = r and 
their number depends on the sign of the discriminant 



D L = {R 1 /\t\ 2 + l) 2 + a\-l. 



(100) 



For Dl < 0, Nq = 1 and the point is of the lemon type. In the opposite case with 
D L > 0, the C points can either be stars or monstars depending on the index, I c = —1/2 
for stars and Ic = 1/2 for monstars. 

As it can be seen from Eq. (91), the derivative Pi is positive if Ic = — 1/2 and, as a 



result, Dl > 0. So, in agreement with the known result, our criterion leads to the conclusion 
that the C points with Ic = —1/2 are always stars. 

At Ic = 1/2, Pi is negative and the point will be monstar provided the gradients |Pi| 
and | ai | are sufficiently large. For example, the points where Pi < — 2\t\ 2 are monstars. As 



is illustrated in Fig. |2(a) in sufficiently thick cells, the C points with k > 2 and Ic = 1/2 
are monstars, whereas the C points in the circle of the smallest radius (k = 0) are lemons. 

More generally, monstars will be the dominating type in regions of pronounced depen- 
dence of the transmission coefficients on the incidence angle. By contrast, a mechanism 
suppressing strong inhomogeneity may drastically increase the fraction of lemons. For ex- 
ample, in liquid crystals, strong director deformations (and large gradients) are inhibited 
due to the elastic energy costs and, thus, disclinations of the monstar type are generally 
unstable pp. 
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For the homeotropic configuration, the condition (84) can be simplified giving the equa- 
tion 



lm(t x t*)sm2ip = 



(101) 



that describes the L lines. From Eq. (101 ) there are two straight lines of linear polarization: 
= ip p and = if) p + 7r/2, where the polarization vectors of incident and transmitted waves 
are parallel, P = ip p . 

Other L lines are circles separating the circles of C points. The radii of the circles, p~^\ 
can be found as the solutions to the equation 5 = nk with k — 1, . . . N — 1. When k is 
even and \t x \ ~ \t y \, from Eq. (93) it can be concluded that, similar to the straight L lines, 



(f) p ip p . If k is odd, the polarization vector rotates with the azimuthal angle of the incidence 
plane and 



ipp + 2(f). These can also be seen in Fig. 2 (a 



C. Director tilt and biaxiality induced effects 



The analytical results presented in the previous section completely characterize the trans- 
mission angular pattern of the light polarization for the homeotropic orientational structure 
in terms of the polarization singularities. In Fig. 2(a), this polarization resolved cono- 
scopic pattern is depicted as the field of polarization ellipses in the projection plane where 
x = tan 9 inc cos inc and y = tan 9 inc sin inc . 

7r/4, the C points are 



Referring to Fig. 2(a) [see also Eq. (89) 



when \t x \ ~ \ty\ and a 



arranged in chains formed by four rays along which they alternate in sign of the handedness 
and of the index. The symmetry of the chain structure can be deduced from the symmetry 



relations (90). 



In particular, the structure is invariant with respect to the operation of reflection: ip — > 
—ip and v — > —v. It additionally posses the symmetry center located at the origin and 
associated with the transformation of inversion: — > + tt. Clearly, each of the two 
straight L lines, <p = if) p and = ipp + 7r /2, is the symmetry axis and the symmetry center 
is also the center of the L circles. 

Algebraically, the symmetry properties of the polarization pattern follow because the 



transmission matrix for the homeotropic cell (67) is diagonal and the transmission coeffi- 
cients, t x and t y , are independent of the azimuthal angle of the incidence plane, = 0i nc . 
For tilted director configurations with < 6& < vr/2, this is no longer the case. 



In this case, the transmission matrix (61) is non-diagonal and its elements depend on 
0d — 0. Generally, there are two types of the azimuthal angle dependent terms: (a) the terms 
proportional to cos2(0d — 0) and (b) the inversion symmetry breaking terms proportional 
to either sin(0 d — 0) or cos(0d — 0). 

In Fig. [2j we show what happen when the tilt angle 9 a varies from #d = to 9d = tt/2, so 
that the orientational structure changes from the homeotropic configuration to the planar 



one. From Fig. 2(b) it is seen that, at small incidence angles, the symmetry breaking terms 
result in the shift of the symmetry center. When the incidence angle #d increases further, 
the center of symmetry eventually leaves the region specified by the aperture leading to the 



single-chain structure of the C points shown in Fig. 2(c 

In the case of the planar structure with d z = (6> d = vr/2), it is not difficult to see 
that the terms dependent on the azimuthal angle are proportional either to dy or to d\. 
So, similar to the homeotropic cell, the origin is the symmetry center for the polarization 
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pattern of light passed through the planar cell. Fig. 2(d) demonstrates that, by contrast to 
the homeotropic structure, where the polarization field is of the elliptic type (see Fig. 2(a)), 
the angular pattern of polarization is of the hyperbolic type for the planar structure. 

The results computed for weakly biaxial orientational configurations are presented in 
Fig. [3| In these calculations, the Euler angles, 0d and 7d, that specify orientation of the 



principal axes ^ are taken to be zero and 7r/4, respectively. In Fig. 3(a), similar to the 



homeotropic cell, the tilt angle 9d is zero, whereas the polarization pattern for the tilted 
biaxial structure [see Fig. |3(b)| w as calculated at 6>d = 10 degrees. From the patterns 
plotted in Fig. |3(a)| and Fig. 3(b) it is clear that biaxiality induces additional deformations 



of the polarization ellipse field as compared to the case of uniaxially anisotropic cells. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the polarization resolved angular patterns of light trans- 
mitted through uniformly anisotropic NLC cells. Such patterns can be considered as the 
polarization structure underlying the conoscopic images measured in experiments with two 
crossed polarizers. 

Experimentally, this structure can be investigated using a suitably modified method of 
the Stokes polarimetry designed for accurate measurements of the Stokes parameters (typi- 
cal experimental setups are described, e.g., in Refs. [251 ED]). The predictions of our theory 
for the homeotropic cell was successfully compared with the polarization patterns experi- 
mentally obtained by the group of experimentalists at Institute of Physics of NASU (Kyiv, 
Ukraine). These results are the subject of a separate joint publication. 

The exact solution to the transmission (reflection) problem for the uniformly anisotropic 
cell was derived in Sec. [IT] by using the 4x4 matrix formalism in combination with the 
orthogonality relations. We deduced the relation linking the operator of evolution and the 



transmission (reflection) matrix [see Eq. (24) and Eq. (26)]. Then the basis of eigenwaves was 



used to construct the operator of evolution (57). After simplifying algebraic rearrangements 



made with the help of the orthogonality relations (18), the final result [see Eq. (61) and 



Eq. (62)] was expressed in terms of 2 x 2 matrices. 

For the case of one interfacial (reflecting) surface, the explicit formulas in the 2x2 matrix 
form were previously derived in Ref. [52l [53] . Algorithms and useful approximations for 
inhomogeneous z-stratified anisotropic media were suggested in more recent papers [Ml ESI 
I5"B"1 15T] . For our purposes, however, it is suffice to have the exact expressions in the compact 
convenient form that emphasize the symmetry properties of the transmission matrix and 
enormously simplifies numerical calculations. An extended discussion of generalizations and 
approximation schemes is beyond the scope of this paper. 

The analytical results were applied to characterize the polarization resolved angular pat- 
terns in terms of the polarization singularities such as C points and L lines. The case of the 
homeotropic director configuration was analyzed in detail. We found that the singularities 
are symmetrically arranged and deduced the simple formula for the index of C points. We 
have also formulated the quantitative criterion to determine the type of the C points. 

The C points are shown to form chains along which they alternate in sign of the hand- 
edness and the index, Iq. At Ic = — 1/2, all the C points are stars, whereas the C points 
with Ic = +1/2 can either be lemons or monstars depending on how strong the depen- 
dence of the transmission coefficients on the incidence angle is. In sufficiently thick NLC 
layers (d > 10 /im for 5CB), the C points of the lemon type are found to be located in the 
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immediate vicinity of the origin. 

The patterns for the tilted, planar and weakly biaxial configurations were computed 
numerically and presented as the polarization ellipse fields. The primary goal of these calcu- 
lations is to demonstrate the effects of director orientation and biaxiality on the polarization 
resolved angular patterns. At this stage we can only hope that these effects can be used to 
develop an improved experimental technique for detecting the peculiarities of liquid crystal 
orientational structures. But our results can be regarded as a first step in this direction. 

Applicability of our results is not restricted to the case of nematic cells that were discussed 
as technoligically important materials. In our considerations, a nematic liquid crystal can 
be easily changed for another uniformly anisotropic medium. 

Our concluding remark concerns the effects related to the polarization state of the incident 



light. For definiteness, in Sec. Ill B and Sec. Ill C this wave was assumed to be linearly 



polarized. More generally, ellipticity and polarization azimuth of the incident wave can 
be regarded as additional governing parameters that may have a profound effect on the 
polarization angular patterns. For example, if the incident light is circular polarized, the 
angular pattern for the hometropic cell is characterized by the only C point located at the 
origin. 

The effects of the incident wave polarization can be studied in terms of generalized po- 
larization patterns in differently parameterized planes. We shall extend on this subject in 
our subsequent publications. 
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